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^Nj ' Abstract. We obtain estimates on the continuous dependence on the coef- 

ficient for second order non-hnear degenerate Neumann type boundary value 
problems. Our results extend previous work of Cockburn et.al., Jakobsen- 
Karlsen, and Gripenberg to problems with more general boundary conditions 
and domains. A new feature here is that we account for the dependence on 
the boundary conditions. As one application of our continuous dependence 
results, we derive for the first time the rate of convergence for the vanishing 
viscosity method for such problems. We also derive new explicit continuous 

P^ ' dependence on the coefficients results for problems involving Bellman-Isaacs 

^A ' equations and certain quasilinear equation. 

^' 

1. Introduction 

In this paper we will derive continuous dependence estimates for the following 
boundary value problem: 

rn ; F{x,u,Du,D^u) = in n (f7cM^), (1.1) 

G{x,Du)^0 on dfl, (1.2) 



m 

where u is the scalar unknown function, Du and D^u denote its the gradient and 
C"^ ' Hessian, and H is a bounded, smooth [W^'°°) domain in M^. Informally speaking, 

^^ , by continuous dependence estimates we mean estimates of the type 

O: |kl-U2|l<|li^l-i^2|| + ||Gi-G2|| 

where ui and U2 are solutions of two different boundary value problem with data 
K^ , Fi,Gi and F2,G2- The exact statement is given in Section [21 

$H ' Equation (jl.ip is degenerate elliptic, (possibly) non-linear, and increasing in u. 

^ This means that the possibly non-linear function F{x,r,p,X) satisfies 

F{x,r,p,X) <F{x,s,p,Y) for aU r < s, X > Y, 

where xSfi, r, sGM,pG M^, and X,Y ^ §^. Here §^ is the set of real symmetric 
N X N matrices and X > in S^ means that X is positive semi-definite. The 
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boundary condition (|1.2p satisfies the Neumann type condition that G is strictly 
increasing in p in the normal direction: 

G(x,p + tn{x)) > G{x,p) + vt, 

for some u > Q and all t > 0, x G d^, p, and outward unit normal vectors n(x) to 
dft at X. The other assumptions on F and G will be specified later. 

The class of boundary conditions G we treat in this paper includes the classical 
Neumann condition, ^ = g{x) in dil, oblique derivative conditions, and non-linear 
boundary conditions like the capillary condition 

^^e{x)(i + \Du\Y^^ in an, 

on 
and the controlled reflection condition 

sup {ja{x) • Du — ga{x)} — in dfl. 

aeA 

In this paper we will assume that g, 9, ga, 7q are Lipschitz continuous functions, 
that 1^1 < (^ < 1 and "/a ■ n > u > 0, and that A is a compact metric space. 

The main class of equations that our framework can handle are equations sat- 
isfying assumption (H2) in the next section. Loosely speaking, this is the class of 
equations where the non-linearity F{x,r,p, X) is uniformly continuous in r,p,X 
locally uniformly in x. This case will be referred to as the "standard case" in the 
rest of this paper. Assumption (H2) excludes most of quasilinear equations, but 
contains fully-nonlinear equations like the Bellman-Isaacs equations from optimal 
stochastic control and stochastic differential games theory: 

inf snp{~tr[a'^'^'''{x)D^u]-b'''^'^'{x)Du~c'^'^'^'{x)u-f''^'{x)}=0 in n, 

where 0i, 02 C M™ are compact metric spaces, c^^^^^ > A > 0, the matrices a^'-'^^ > 
0, and the coefficients are Lipschitz continuous uniformly in ^i, 02- In Sect ions [3] and 
[SI we give all the details, more examples, and extensions to problems on unbounded 
domains, time-dependent problems, and certain quasilinear equations like e.g. 

_[l-^——)D^+Xu = .fix) m n. 



-tr 



Since these equations may be degenerate and non-linear, their solutions will in 
general not be smooth. In this paper we work with a concept of weak solutions 
called viscosity solutions, a precise definition is given at the end of this introduc- 
tion. Viscosity solutions are at least continuous in the interior of O. The boundary 
conditions will be interpreted in the weak viscosity sense which essentially means 
that either the boundary condition or the equation has to hold on the boundary. 
This allow us to have well-posed problems even when the boundary conditions 
are classically incompatible. The solutions can realized by the vanishing viscos- 
ity method, and they will be discontinuous at parts of the boundary where the 
boundary conditions are classically incompatible. 

An overview of the viscosity solution theory, including Neumann boundary value 
problems, can be found in the User's Guide [12]. The viscosity solution theory for 
Neumann type boundary value problems was initiated by Lions |29j in 1985 for first 
order equations, and has been developed by many authors since, see [30l [HI [TJ [2j 
[201 [2] and references therein for various aspects of this theory. Today there are two 
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leading approaches, one due to Ishii |18, and another one due to Barles [llE]- They 
apply under slightly different assumptions and will be discussed below. 

Starting with the standard case, i.e. nonlinear equations (jl.ip satisfying (H2), we 
prove under natural and standard assumptions, that these boundary value problems 
have unique Holder continuous viscosity solutions. The Holder regularity results 
are new and extend the Lipschitz regularity result of Barles [1], and we give for 
the first time a complete proof of such a regularity result. We note that these 
regularity results are global up to the boundary. Local up to the boundary Holder 
estimates have previously been obtained by Barles-Da Lio |3] for a different class 
of equations. Whereas our equations are degenererate but strictly increasing in 
the u argument (assumption (H3) in the next section), their equations are weakly 
non-degenerate satisfying some sort of "strong ellipticity condition" but are not 
necessarily increasing in u. The arguments needed to prove the two types of results 
are also different, except for some ideas on the construction of test functions that 
are needed in some of the proofs. 

Next we prove continuous dependence results comparing Holder continuous (sub 
and super) solutions of different boundary value problems. The results we ob- 
tain include both continuous dependence on non-linearities for the equation and 
the boundary condition. The results concerning the dependence on the boundary 
condition are completely new, at least in a viscosity solution setting, while the re- 
sults we obtain for the equations apply to much more general boundary conditions 
(including non-linear ones) than earlier results. 

Continuous dependence results for the type of equations we consider in this paper 
have previously been obtained by e.g. Cockburn et.al. TT], Jakobsen and Karlsen 
[3TJ [55i i23^ , and Gripenberg [17^. In all these papers viscosity solutions methods are 
used. In some cases such results can also be obtained from probabilistic arguments, 
see e.g. [14] for results for Bellman equations set in R^. Papers [211 1221 [23] 
treat very general classes of equations set in R^ or R^ x [0,T), [TT] treats zero- 
Neumann boundary value problem for x-independent equations, and |17j treats a 
zero-Dirichlet boundary value problem. 

In the two last papers the domain fJ is convex and possibly unbounded and in 
the last paper further restrictions on the class of equations are needed (because of 
the Dirichlet condition) and the Dirichlet condition is taken in the classical sense. 
All these papers treat more general quasilinear equations than we can treat here, 
e.g. p-Laplace type equations for p > 2. 

The technical explanation for the differences between our continuous dependence 
result and the above mentioned results lays in the choice of test function we use. 
To handle weakly posed Neumann boundary conditions, the idea is to use a test 
function that will never satisfy the boundary condition. The effect will be that the 
equation holds also at the boundary, and that the classical viscosity solution com- 
parison argument can be used (see the following sections). To achieve this the usual 
test function has to be modified and the extent of the modifications depend on how 
smooth and non-linear the Neumann condition is. To handle possibly non-linear 
boundary conditions or Holder continuous solutions in combination with bound- 
ary reflection directions that are only Lipschitz functions in the space variable, it 
seemed that the only available or at least the most optimal test function to use, is 
the one constructed by Barles in [T]|2] . As opposed to the basic test function used 
in the other papers on continuous dependence, the test function of Barles is not 
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symmetric in its arguments {x and y) and therefore it does not have equal x and y 
gradients. We loose a cancellation property in the comparison proof and hence can 
not handle as general gradient dependence in the equations as with the basic test 
function. In this paper we consider the same class of non-singular (!) equations as 
Barles in [U [2] , and this excludes most of the quasilinear equations considered in 
[m [2T| [22l [TT] , including p-Laplace equations for p 7^ 2 (see also remark ()5.ip in 
section 5.). 

At this point we mentioned that a different test function has been constructed by 
Ishii in [T5]. Compared with Barles, Ishii is able to treat less regular domains but 
with more regular (and less non-linear) boundary conditions (e.g. C^ domains and 
W'^''°° reflections), see [1] for a more detailed comparison. Using Ishii's test function, 
continuous dependence results could probably be obtained under a different set of 
assumptions (see above). We have not considered this case. 

We also point out that we can handle u-depending boundary conditions only 
through additional arguments involving transformations. This is in contrast to the 
general comparison results obtained by Barles p] under similar assumptions for F 
and G. In [2] u-dependence is handled directly by a sort of localization argument 
(Lemma 5.2 in [?) 

which only works when you send some parameter of the test function to zero. In 
our continuous dependence arguments, we will have to optimize with respect to this 
parameter and the optimal choice will in general not be zero or even small. See the 
treatment of parameter e at end of the proof of Theorem 12.21 One way to handle 
M-depending Neumann type boundary value problems, is to transform them into 
problems with no u-dependence, then to use our results, and finally to transform 
back. 

We do not consider such transformations in this paper, instead we refer to [3] 
where such transformations have been considered in a rather general setting. 

Continuous dependence results have to do with well-posedness of the equation. 
Typically the boundary value problem we consider model some physical process, 
and the data is measured data. A continuous dependence result then implies that 
small measurement errors only produce small errors in the solutions. Any reason- 
able model should satisfy such a requirement in particular in view of numerical 
computations. Moreover, continuous dependence results have been used in many 
other contexts. They play a key part in the shaking of coefficients approach of 
Krylov to obtain error estimates for approximation schemes for Bellman equations 
PSI [^ [771 m ini IS] , in Bourgoing 'W and in '^V\ they are used to obtain regularity 
results, and they have been used to estimate diffusion matrix projection errors [7], 
source term splitting errors |24;, and errors coming from the truncation of Levy 
measures [23] . They have also been used to derive the rate of convergence for the 
vanishing viscosity method [HI [211 [221 [17] , see also e.g. [TO] . 

The paper is organized as follows: In the next section we state the assumptions 
on the boundary value problem (|1.1[) and (jl.2p in 

the standard case and give well-posedness and Holder regularity results. We state 
the main result, the continuous dependence result, and as an immediate corollary 
we derive an estimate on the rate of convergence for the vanishing viscosity method. 
The proofs of the main result along with the regularity result are proven in Section 
[3l and in Section [4| we apply our main result to obtain new continuous dependence 
results for boundary value problems involving Bellman-Isaacs equations. We give 



CONTINUOUS DEPENDENCE 5 

several extensions of our results in Section [5l to time-depending equations, equa- 
tions set on unbounded domains, and certain quasilinear equations. Finally, in the 
Appendix we derive the test function used in the proofs in Section [3] along with its 
properties. 

Notations. We let | • | denote the Euclidean norm both in R™ (vectors) and R™^p 
(matrices) for to,p G N. We denote by S^ the space of symmetric N x N matrices, 
tr and '^ denote trace and transpose of matrices, and < denote the natural orderings 
of both numbers and square matrices. For a, 6 G R we define a V 5 = max(a, b) and 
aAb = min(a, b). We will also denote various constants by K or C, and their values 
may change from line to line. 

Let BUSC{U), BLSC{U), C{U), and Wp^°^{U) denote the spaces of bounded 
upper and lower semicontinuous functions, continuous functions, and functions with 
p essentially bounded derivatives, all functions defined on U . If / : M^ -^ ]]j"ixp jg 
a function and a S (0, 1], then define the following (semi) norms : 

|/|o = sup|/(x)|, [/]„^ sup '^f^j{y , and |/U = |/|o + [/].■ 

By C°^"(rJ) we denote the set of functions / :— > M with finite norm \f\a- 
We end this section by recalling the definition of a viscosity solution: 

Definition 1.1. An upper semicontinuous function w is a viscosity subsolution of 
(|l.ip and (|1.2p if for all G C^(f2), at each maximum point xq ^ Ct oi u — cj), 

F{xo,u{xo), D(l){xo),D^(l){xo)) < if xo G ^i, (1.3) 

mm{F{xo,u{xo), D(j){xo), D'^(j}{xo)), G{xo, Du{xq)) < ii x^ e dVL (1.4) 

An lower-semicontinuous function u is a viscosity supersolution of (jl.ip and ()1.2p if 
for all (f) e C^(fi), at each minimum point xq d fl oi u — (f>, 

F{xo,u{xo), D(j)ixo),D^(f>{xo)) > iixoefl, (1.5) 

max{F{xo,u{xo), D(j){xo),D^(j){xo)), G{xo, Du{xo)) > ii xq e dQ (1.6) 

Finally u is a solution when it is both a super and a sub-solution. 

2. The main results 

In this section we consider the standard case (when assumption (112) below 
holds) . Following [1] [2] we state the assumptions on the boundary value problem 
(|l.ip and (jl.2p and give results on comparison, uniqueness, and existence of solu- 
tions. Then we give new Holder regularity results extending the Lipschitz regularity 
result of [1] in two ways: we allow Holder continuous data and small A (see assump- 
tion (H3) below). We also give a complete proof. The main result of this paper, 
the continuous dependence result, is then stated, and as an immediate consequence 
we derive an explicit rate for the convergence of the vanishing viscosity method. 

Here is a list of the assumptions we will use, starting by the domain: 

(HO) il is a bounded domain in R^ with a W^'°° boundary. 
For the equation we use the following standard assumptions: 

(HI) F e C{n X R X R^ X §^). 
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(H2) There exists a modulus lur^k (a continuous, non-decreasing function 
satisfying ujr^k{0) = 0) such that 

F{y,r,q,Y) - F{x,r,p, X) < u;ii^K{\x - y\ + ^\x - y\^ + 7f + e^ + B^ 

for £,77 e (0,1], B > 0, x,y e ^, r £ R, \r\ < R, p,q e M^ and 
X,Y eS^ satisfying \x-y\< Kr]s, \p-q\< K{rf + e^ + B), \p\ + \q\ < 
K{^+rf + £^ +B), and 

X \ ^ K f Id -Id\ ^^,2 2 r^^fld 0\ ,^_,, 
,0 ^Y)^7^[-Id /rfj+^(^ +^ +^)(,0 Id)- (2-1) 

(H2) There exists a e (0, 1] and Kr>0 such that 

F(y, r, g, r) - F{x, r,p, X) < Kr{\x - y]" + 1|x - y\^ + rf + e^ + B 

where e, 77, _B, i?, x, ?;,p, q, X, Y are as in (H2). 
(H3) For every x,p, X, and for any _R > 0, there is A/? > such that 

F{x,r,p,X)- F{x,s,p,X)>\R{r - s) for -R<s<r<R. 
The possibly fully nonlinear Neumann type boundary condition satisfies: 

(HBl) There exists v>Q such that for all /x > 0, a; £ dVl,p e R^, 

G(x,p + ^n{x)) — G{x,p) > v^, 

where n{x) is the unit outward normal at x. 

(HB2) There exists a constant K such that for all x,y E dft and all p, g G 1 

\G{x,p) - G(y, q)\<K[{l + \p\ + \q\)\x -y\ + \p-q\]. 

Remark 2.1. In general there is a trade off between the regularity of the boundary 
d^ and the generality and smoothness of the boundary condition G, see [T] for a 
discussion. (HO) compensates for very general non-smooth boundary conditions. 

Remark 2.2. Assumption (H2) plays the same role as (3.14) in the Users' Guide 
[T^ . By this assumption the equation is degenerate elliptic. Moreover, it is a re- 
fined version of assumption (H5-1) in [2] containing also a new parameter B. In 
the proofs, this parameter will be used to carry information from the boundary 
conditions (which are never satisfied, see the introduction) over to the equations. 
Assumption (H2) is a strengthening of hypothesis (112) which yields Holder regu- 
larity results. 

By (H3) the equation is strictly increasing in the u argument. Assumption 
(HBl) is the Neumann assumption, saying that the boundary condition G, contains 
non-vanishing and non-tangential (to 9f2) oblique derivatives and it is a natural 
condition to insure the well-posedness of the problem. 

We now state a comparison, uniqueness, existence, and regularity result for so- 
lutions of (fTTTI) and m 



)A' 



Theorem 2.1. // (HO), (HI), (H2), (H3), (HBl), and (HB2) hold, then the fol- 
lowing statements are true: 

(a) If u is a BUSC{n) subsolution and v is a BLSC{fl) supersolution of (jl.ip and 
(11.21). then u < v in fl. 



CONTINUOUS DEPENDENCE 7 

(h) If \r in (H3) is independent of R, then there exists a unique solution u G C{Vt) 
of dTU) and (fO]) . 

(c) Assume (H2) also holds, u G C(n) is the solution of (jl.ip and (jl.2p . and 
A :— \u\o > 0- ^/len there are constants (3 G (0, a] anrf X (^onZy depending on the 
data and X) such that 

\u{x) — u{y)\ < K\x — yf in Q x Cl. 

Furthermore, there exists a constant A > (only depending on the data) such that 
if X > X then (3 — a (the maximal regularity is attained). 

The comparison principle in (a) correspond to Theorem 2.1 in [2]. The unique- 
ness part in (b) follow from (a), and existence follows from Perrons method [19] 
since w{x) := M — Kd{x) is a supersolution of (jl.ip and —w is a subsolution of 
(|1.2p . if Af, K >{) are big enough, and d is the W^'"^ extension of the distance func- 
tion defined in the Appendix, see Section 4 in [^ for similar results. The regularity 
result, part (c), will be proved in Section [31 

Remark 2.3. The regularity results in part c) are global up to the boundary. Local 
up to the boundary Holder estimates have been obtained by Barles-Da Lio |3] using 
different techniques and assumptions on the nonlinearity of the equation. See the 
introduction for a discussion. 

Now we proceed to the continuous dependence result. We will derive an upper 
bound on the difference between a viscosity subsolution ui of 

Fi{x,ui{x),Dui{x),D'^ui{x)) = in f}, (2.2) 

Gi(x, Dui{x)) 
and a viscosity supersolution U2 of 

F2{x,U2{x),Du2{x),D^U2{x)) = in 0, (2.3) 

G2{x,Du2{x)) 

We assume the following estimates on the differences of the two equations and of 
the two boundary conditions. 

(Dl) There are 6i, §2 > 0, and Kf{K) > such that for any K > 0, 

F2{y,r,q,Y) - Fi{x,r,p,X) < Kp{K)[f]^ +di + ^dl + By 

for < £ < 77 := e^^ < 1 with a ^ a A l3, B > 0, x,y G Cl, r e R, \r\ < K, 
p,qeR^ and X,Y gS^ satisfying |a; - y| < K-qe, \p-q\< Krf^ + KB, \p\ + \q\ < 
K{^+r]'^ + B), and 

X \ ^ K f Id -Id\ ^^, 2 ^^ (Id 

(D2) There are ^i, ^^2, ^g > such that for all a; G 9^7 and p e M^, 

G2(x,p) - Gx{x,p) < Kcini + ^l2\p\)■ 

Remark 2.4. Assumption (Dl) is a "continuous dependence" version of (H2) and 
(H2) in this paper, and assumption (3.14) in the Users' Guide [12]. A similar 
assumption is used in Theorem 2.1 in [22j . 

By P and a we denote the Holder exponents of the solutions and data respec- 
tively. In general a > P, and equality only holds when A in (H3) is big enough. 
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Since \x — y\ < Kei], r/ ~ e^^ imply '^ i < Krf' and \x — y\^ < Kjf, the 
Fi ~F2 inequality in (Dl) will be implied by the following more standard inequality 

F2{y, r, g, Y) - Fi{x, r,p, X) < k{\x - y|" + l|a; - y\^ + Si + ^S^ + ^^ + b). 

Remark 2.5. on assumption (D2). In the case of oblique derivative boundary con- 
ditions, Gi{x,p) — ji{x) ■ p — gi{x), i = 1,2, and 

\G2{x,p) - Gi{x,p)\ < \{gi - g2) + |o + l7i " 72|ob|. 

Our main result is stated in the following theorem: 

Theorem 2.2 (Continuous Dependence Estimate). Assume (HO), (HI), (H3), 
(HBl), and (HB2) hold for Hi,H2,Gi,G2, and ui,U2 £ C^^^in) for (3 £ (0,1]. 
Define v^ = (y^ V V2){vi A 1/2) and X = Aij^^i^ V A2,|„2io. 

// (Dl) and (D2) hold and ui and U2 satisfy the boundary value problems (12. 2|) 
and (j2.3p respectively, then there exist a constant C > (depending only on Kp, 
K, Kg , \ui\f3, \u2\13, a., (3) Such that 

Amax(ui--U2) < Cdi+dn H h — 

n V V \ V J J 

Remark 2.6. As far as we know this is the first result giving continuous dependence 
on the boundary condition. The result also extends the earlier continuous depen- 
dence on the equation type of results of [TTJ [^TJ [551 EZ] since much more general 
boundary conditions arc considered (but at the expense of less general equations!). 

We prove Theorem 12.21 in Section [S] An immediate consequence of this result 
is an estimate on the rate of convergence for the vanishing viscosity method. For 
/U > we consider the solution u^^ of 

F{x,u,Du,D'^u) = fiAu in O, (2.4) 

with bomidary condition (II. 2p . The result is the following: 



Theorem 2.3. Assume (HO), (HI), (H2), (H3), (HBl), (HB2), p > 0, and that 
u and u^ solve (jl.ip / (|1.2p and (|2.4p / (|1.2p respectively. Then u and u^ belong to 
C^'f^in) for some f3 e (0,a] and 

\u-u^\o<Cp''/^ 

Proof. Regularity follows from Theorem 12. II By assumption (H2) 

[F{y, r, q, Y) - ptiY] - Fix, r,p, X) <C{\x- y^ + ^^ + if + e') - pirY, 

and inequality (|2.ip implies that — trF < C-^ + small terms. Theorem 12.21 imme- 
diately gives u — U/^i < C^i^''^. A lower bound can be found in a similar way. D 

Remark 2.7. This result seems to be the first such result for complicated boundary 
condition. We refer to [HI [17] for results on weak 0-Neumann or classical 0-Dirichlet 
problems, to [31] for results on linear Neumann boundary value problems for first 
order equations, and to [HI [22] for result in R^ or (0,T) x M^. 



Remark 2.8. The vanishing viscosity method has been studied by many authors 
dealing with weak solutions of nonlinear PDEs. The method has been used to 
obtain existence (and uniqueness!) of solutions for degenerate (e.g. first order) 
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problems by taking the limit a,s fi —> (see e.g. [9, 32 J. and it is well-known that 
it is strongly related to the problem of proving convergence rates for numerical 
approximations of such problems (see e.g. [T51[ST| ). 



3. Proofs of Theorems 12.21 and 12. II (c) 
Proof of Theorem \2.2[ First we assume without loss of generality that 

0l,02, , < 1. 

V V 

If this is not the case then the theorem holds since 

u,~u2< (luilo + \u2VM5i + <5f + ^ + f ^ 



where a = a /\ (3. Then we double the variables and consider 

'il}{x,y) ^ui{x) -U2{v) ~ (j){x,y) and M ^ u\a.yi_ip{x,y) = ip{x,y), 

where for A, _S > 0, 

(/.(x, y) ^ l^\x - y\^ + ^ (d{x) - d{y)f - B{d{x) + %)) 

^ ,x + y 2{x~y) 
2 e 

and C2{x,p) = C2,a{x,p) with a = ije = e^^^ [rj = e^^^ by (Dl)). The functions 
6*2,0 and d are defined in the Appendix, and the smooth function (j) was introduced 
by Barles in [2J . We refer to the Appendix for the proofs of the properties of (p. 

The existence of a point (x, y) follows from compactness of 51 and the continuity 
of all functions involved. Since {x,y) is a maximum point, 

2'0(S, y) > ^{x, x) + ilj{y, y). 
Moreover, if A is big enough. Lemma [A. 31 of the Appendix implies that 

0(5;, y)>^\x- y\^ - K^e^ - B{d{x) + d{y)), (3.1) 

and Holder regularity of ui and U2 combined with the last two inequalities yield 

^\x - y\^ < K^\x - yr y e^ 

for some constant Ki depending on Kq and the Holder constants of ui and U2 (but 
not on B). Equivalently, since 77 — e^^^ by (Dl) and e < rj, 

\x - y\ < kie^ = ki-qe and -^\x ~ y\^ < Kie^ ^ KiT]'^ . (3.2) 

Now we choose A and B in the test function (p to insure that when x ov y belong 
to the boundary 951, then the boundary conditions can not hold there. See Lemma 
IA.4I of the Appendix. This means that the equations always has to hold at x and 
y. The precise choices of A and B are 

B = K{7f + s'^) + —(fii+^i2-] and A = K, 
for some K only depending on the data of the problem. 
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By the maximum principle for semicontinuous functions, Theorem 3.2 of the 
sers' guide" 12j, there are {p,X) e J^^ui{x) and {q,Y) £ J^~U2{y) such that 

p = D^(f>{x,y), q = -Dy(t){x,y), 

Q ^y] < [Id + e^D^^ix, y)]D^^{x, y). 

Using the definition of viscosity sub and super solutions at x and y (and Lemma 
IA.4|) we get 

Fi{x,ui[x),p,X) <Q < F2{y,U2{y).q,Y). 
We rewrite this as 

Fi{x,ui{x),p,X) - Fi{x,u2{y),p,X) < F2{y,U2{y),q,Y)) - Fi{x,U2{y),p,X). 

(3.3) 
By Lemma [A. 5| the definitions of p, q, X, Y, and £ < ?7 < 1, it follows that 

\p-q\< Krf + 25, 
X Q \ ^ K f Id -Id\ ^^, 2 n^ [Id 

again for some K only depending on the data of the problem. Since we also have 
(|3.2p . we are in a position to use assumption (Dl). So if ui{x) — U2{y) > 0, then 
(Dl) and (H3) applied to (gS]) yield 

\i{ui{x) - U2m < Kf{K) (r/2 + 5i + I52 + S) . 

By (|3.ip and the definition of xjj, it follows that 

ui{x) ~U2{x) < il^e{x,x) < ips{x,y) < ui{x) -U2{y) + 2B{d{x) + d{y)). 
Therefore the two previous inequalities and the choice of B implies that 

Ai(ui(x) - U2{x)) < K(f + 5, + l^i + 1^^1±^) . 

Remember that rj — rj{e) — e^^^ and let £1 and £2 be defined by 

V{eif = \5l or v{eif^5l 

,, ')(e2) g 

?7 £2 ' = 2^ or 7/ £2 2 = ^. 

Now with £ = £1 V £2 (< 1 by assumption) it follows that 

\i{ux{x) ~ U2{x)) <K(^Si+6'i + ^ + ^ 

A closer look at the proof reveals that we may replace Ai by Ai V A2. □ 

Proof of Theorem \2.1\ (c). We start by proving a-Holder regularity when A is big 
(the last statement of Theorem l2.1l (c) ) . The proof is similar to the proof of Theorem 



CONTINUOUS DEPENDENCE 11 

12.21 except that we have to modify the test function and use a bootstrap argument. 
The modified test function is 

_c.(i±i. ^--''-"'";:'"''-'') )K.)-.(,)) 

- Keia + e^){d{x) + d{y)). 

We refer to the Appendix for the definitions of Ca and d. Playing with the parameter 
a, we will use a bootstrap argument to prove that u has the right regularity. 

The new test function satisfies similar estimates as the ones given in Lemmas 
IA.3I - IA.5I The moral is that the new terms coming from the exponential term are 
not worse than the old terms. We refer to p for such estimates given in the full 
generality (but with a different choice of a). 

Now let £ < 1 and double the variables defining 

M -.^ -tpix^y) = sup^p{x,y) where 'ip{x,y) ^ u{x) - u{y) ~ (j)a{x,y). 
If A is big enough, (an easy extension of) Lemma fA.SI and the inequality 2ip{x, y) > 



(3.4) 



tp{x, 


x) + 


ip{y, y) imply 


thai 










1 
2e2 


e-^=' 


{d(x)+d{y 


))|x- 


-V? 


< 2[uix) 


- u{y)] + Koe^ < 


Define ry 


= K- 


-^X\x- 


y? 


with K = e^^' 


{2\u\o + Ko). By 










ri' 


< 1 


and \x 


-y| < K^^^v^- 



We proceed as in the proof of Theorem [ 

By arguments similar to the ones in the proof of Lemma IA.41 if A, K^, and 
Kb are big enough (not depending on e, a or B), then the equation holds even if 
{x,y) lies on d{il x fi). Compared with the proof of Theorem 12. 2[ the exponential 
allows us to cancel at the boundary all terms of the form ■ir\x — y\'^ = Ki]^ and use 
B = Ksia + £2^^) at each step. 

Note that D(j)a and D'^cfia still satisfy inequalities (|A.12p and (JA.16P in Lemma 
IA.5I We will choose a such that inequality (jA.16|l takes the form of (|2.ip . i.e. we 
choose a such that -rj'^ < K. Since 77^ < 1 we choose a — e. Again we use the 
definition of viscosity solutions and subtract the equations (inequalities) at x and 
y using the maximum principle for semi continuous functions. By the appropriate 
version of Lemma I A. 5 1 and the definition of rj^ and B we can now use (HI) and 
(H2) to get 

Xiuix) - u{y)) <k{\x- yl" + -^\x - y\'' + v^ + e^ + Keie + £5^)) . 

By Young's inequality, the definition of i]^ , and £ < 1 we have 

X{u{x) ~ um < k[^\x - yf + e'^'^y 

When A is big enough, an appropriate version of Lemma [A. 31 the definition of M, 
and < d < 1, imply that 

u{x)-u{y) = M+(f>ix,y) > M+—e-^^''\x-y\^-KoS^-KBie+e^){d{x)+d{y)). 
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Combining the two last inequalities and using that e < 1 leads to 
\M <(K~ L-2K.)1\^ _ y|2 + Ke^^^ 

2 £^ 

If A is big enough, XM < Ke^^^^" , and the definition of M leads to 
u{x) - u{y) - Mx,y) <M< ye'^*^ 

for every x,y G il. Now by the definition of (pa, the properties of the distance 
function, and Young's inequality, we have 

1 



u{x) - u{y) < K-^\x - yf + Ke^'^^- 



If |a; — y| < 1 we may take e = |a; — y| 3 when 1 < ^^ and e^-a = |2;_y|" otherwise, 
the result (since we may also interchange x and y) is that 

Hx)-u{y)\<K\x-y\i'^'^. (3.5) 

If \x — y\> 1, the result still holds since then \u{x) — u{y)\ < 2|u|o|a; — y| a'^". We 
are now done if a < | . 

If a G (§,!], we restart the proof using the regularity estimate p.Sp to get a 
better choice of a such that -ry^ < K. From (|3.5p and the first inequality in (13. 4p . 

2 
?7^ = i^^^ — |x - yp < i^|;r - 2/|3 V £^ and hence ri<Ke^'3\/e, 

so the new choice of a should be er/^ — e^ V e^. But this quantity is less than e"^ 
so we may instead take a = e^ which still implies -r/^ < K. Now it is a simple 
exercise to redo the proof and show that for A big, 

\u{x) — u{y)\ < K\x — y^^ for x,j/Gf2, 

and this completes the proof of the last part of Theorem 12.11 

Now we will prove the first part of Theorem 12.11 (c) using the result we proved 
above and an iterative argument of Lions [55]. Here we only sketch parts of the 
argument, since the details can be found in [22] for similar equations. The idea is 
to consider for /i > 

F{x, u"+\Du"+\D^u"+^) + fiu"+^ = fiu" 

with boundary conditions (jl.2p and noting that u" converge uniformly to u. If 
fi is big enough the above proven result applies, and a careful look at the above 
argument reveals that when A + /i > K{= A), then 

\u''+\x) - u''+\y)\ < ( ^ ^*""*" + other terms') \x - yl", x,yen,\x-y\< 1. 
VA + /i — A / 

Furthermore the comparison principle yields 

/i + A \fi + \/ 

When |a; — y| < 1, the rest of the proof is exactly as in ^2] and we omit it. When 
\x — y\ > 1 any Holder estimate holds since u is bounded. The result is a Holder 
estimate for any A > 0, but with a Holder exponent that is smaller than a. D 
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4. Bellman-Isaacs type boundary value problems 

In this section we apply our results in Section [2] to Bellman-Isaacs equations and 
several different types of boundary conditions. The Bellman-Isaacs equations are 
of the form 

inf sup {-tT[icrcr'^f^'^'{x)D^u] - b^^'^'{x)Du - c^^'^'{x)u - f^-^^{x)\ = 

(4.1) 

in f2. Assumptions (HI), (H2), (112), and (H3) are satisfied [T2t[22] if we assume: 

1. (T^i'^2 g^jjjj ipifii j^j-g Lipschitz continuous in x uniformly in 9i,92, 

2. c^i'^^ and f^^-^^ are a-H61der continuous in x uniformly in ^1,6*2, 

3. c^i'^2(x) > A > for all x, 61,62, and 

4. 01, 02 are compact metric spaces. 

Next, we list some typical boundary conditions we can consider: 

(a) The classical Neumann condition: 

7— — g(x) in dfl. 

(b) The oblique derivative condition: 

—- — qix) in do.. 

(c) The capillary boundary condition: 

-2!:^0{^x)(l + \Du\^f/^ in dn with |^(x)| <w<l. (4.2) 

(d) The "controlled" reflection boundary condition: 

inf snp{-i^^^^^{x)-Du-g^^^^^{x)}^{) in dO.. (4.3) 

"£©1/3692 

Here n{x) is the outward unit normal to dVl. Assumptions (HBl) and (HB2) hold 
in all cases if assumption 4 holds along with 

5. there exists v > such that 

7(2;) • n{x) > V and 7^' ^(z) • n(x) > v uniformly in 61,62, 

6. 5, 7, ^, 7^^'^^, g^^'^^ are Lipschitz continuous in x uniformly in 61,62- 

Now we state new continuous dependence results for the for the Bellman-Isaacs 
equations (|4.ip combined with the controlled reflection boundary conditions (14. 3^ : 

Theorem 4.1. Assume ui and U2 satisfy the boundary value problem (|4.ip and 
(14. 3p with coefficients cri,bi,ci, fi,ji,gi and 172, 62, C2, /2, 72,52 respectively, where 
both sets of coefficients satisfy assumptions 1-6 above. 
Then ui,U2 belong to C'^'l^iCl) for some /3 G (0, a], and 

\ I \ ^ /~^ ^ I ^li^2 ^1-^21/3 I \ J 0-1, 62 1,0-1 ,62 \0 

X\ui-'U2\a<C sup ki -^2'- 'lo + |0i ■ ' -O2'' '|5 



61x62 

C sup 
61x62 



C2 lo + 1/1 - 72 lo 



— sup |.9i'''-ff2 I0 + — sup |7l -72 lo- 

V 61X62 ^ 61X62 
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This result is a direct consequence of Theorems 12.11 and 12.21 In this case Si 
correspond to the second hne in the estimate, 

r2 ^ .,„T^ ri^^i-^2 ^^1:^2 12 I \jdi,02 i,'5l,^2|2i 

d2 = O sup [\(J^ -0-2 I + |0i -02 I J, 

61X62 

Ml = sup |5i -52 lo, A*2 = sup |7i -72 lo- 

61x62 61x62 

The dependence on the equation is as in [3TJ [52] and the derivation of 5i and 82 is 
explained there. 

By Theorem l2.3l we have for the first time the rate of convergence of the vanishing 
viscosity method for the boundary value problem ()4.ip and ()4.3p . i.e. 

inf sup {-tT[{(ja^)''^-^^{x)D^u] - h^^-^''{x)Du - c''^-^^x)u - /''^'(a;)} = /lAu 
eieGi 82662 

(4.4) 

in ri, with (|4.3p as boundary conditions. The result is the following: 

Theorem 4.2. Assume u and u^ satisfy (j4.ip and (|4.4p respectively with boundary 
values (j4.3p . and t/iat assumptions 1-6 hold. 

Then u, u^ belong to C^'^{Cl) for some (3 G (0, a] and 

\u-Ui_,\o < C/i2. 

5. Extensions 

It is possible to consider many kinds of extensions of the results in this paper. We 
will consider three cases: (i) fl unbounded, (ii) time dependent problems, and (iii) 
quasilinear equations. In the two first cases the results cover e.g. Bellman-Isaacs 
equations under natural assumptions on the data. 

5.1. Unbounded domains. Let il be unbounded and let (HOu) denote assump- 
tion (HO) without the boundedness assumption. If we assume that our sub and 
supersolutions u and v are bounded, then we will get continuous dependence and 
regularity results simply by following the arguments in this paper replacing the test 
function (f)a by the standard modification 

Mx,y) + ii\x\^ + \y\''), 7>o. 

The new test function will insure existence of maximum points when we double 
the variables, and at the end of the proof it turns out (as usual) that all terms 
depending on 7 will vanish when 7 — > 0. In the proof B will now depend also on 
the 7-terms and the 7-terms will tend to zero as 7 -^ with a speed depending 
on B, see assumption (Dlu) below. By careful computations and fixing e before 
sending 7 — > we can conclude as before. We refer to [22^ for the details when 

The corresponding continuous dependence result will now be given without fur- 
ther proof. We modify assumption (Dl) so it corresponds to our new test function, 
see also |22) : 

(Dlu) There are ^1, (52 > 0, a modulus uj, and Kp{K) > 0, such that for any 
K>0, 

F2iy,r,q,Y) - Fi(x,r,p,X) < KF{K){rj^+Si + ^Sl + B + j{l + \x\^ + \y\^) 
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for £,7 e (0,1], 7] := e^^, B > 0, x,y e H, r e M, \r\ < K, p,q e M^ and 
X,Y £ 3^ satisfying \x ~ y\ < K-qe, \x\ + \y\ < j^/'^uj {-/){! + B), \p - q\ < 
K{r^^ + B + 7(|x| + Iz/D), \p\ + \q\ < if (f + ry2 + B + ^(\x\ + \y\)), and 

X Q \ ^ K f Id -Id\ , ,w 2 . - . ^ fid 



-Y)^ VA-Id wj+^('^ +^ + ^\0 Id 

Theorem 5.1 {VL unbounded). Assume (HOu), (HI), (H3), (HBl), and (HB2) hold 
forHi,H2,Gi,G2, andui,U2 € C°^^{n) for (3 C, (0,1]. Definev^ = {viViy2){i^i/\i^2) 
and X = Xi^\ui\o V A2,|„2|o. 

// (Dlu) and (D2) hold and ui and U2 satisfy the boundary value problems (12. 2p 
and (12. 3p respectively then there exist a constant C > (depending only on Kp, 
K, Kq , \u\\ji, \u2\fi, a) such that 



Amax(ui-U2) < Coi+O, H h — 



n 
5.2. Time dependent case. Consider a Cauchy-Neumann problem of tlie form: 

Ut + F(t,x,u,Du,D^u)^0 in {Q,T) x n, (5.1) 

G{x,Du) = on (O,r)x90, (5.2) 

u{0,x)^uoix) on {0} X 17. (5.3) 

In this case we get results by similar arguments as above by replacing the test 
function (pa by 

at + e'^*'(j)a{x,y), a > Q. 

We have to replace assumptions (HI) - (H3) and (Dl) by assumptions (Hip) - 
(H3p) and (Dip) depending on t. In (Hip) we assume in addition continuity in i, 
in (H3p) we allow A^ > 0, and in the last two assumptions ((H2p) and (Dip)) we 
simply assume that (H2) and (Dl) hold uniformly in t. Note that one can always 
reduce a problem with A/j G M, via an exponential scaling of u, to a problem with 

Xb,>G- 

Now existence, uniqueness, and regularity results follows as before by appropri- 
ately choosing the constants a and K. Note however that in the result correspond- 
ing to Theorem 12.11 (c) the Holder exponent is always a and "maximal regularity" 
is achieved regardless of the value A. We refer to [H] for such results in the case 
fi = MJ^ . Now we state the continuous dependence result without further proof. 

Theorem 5.2 (Time dependent case). Assume (HO), (Hip), (H3p), (HBl), and 
(HB2) holdforHi,H2,Gi,G2, ui,U2 S C([0,T] x Cl), anduifl,U2,o € C^^^'in) for 
some a € (0, 1]. Define v^ = (y\ V V2){y\ A V2)- 

If (Dip) and (D2) hold and Ui and U2 are sub and supersolutions of initial 
boundary value problems (|5.ip . (I5.2[) . and (|5.3[) respectively for Fi,Gi,ui.o and 
F2,G2,U2fi, then there exist a constant C > (depending only on Kp, K, Kq, 
|'"i,oU>l''^2,o|a7 T, a) such that for t £ (0,T), 

max(ui(t,-)-M2(i,-)) < K^i.o - W2,o)+|o + C7t(5i + 5? + ^ + f^)") . 

Note that we do not need to assume that u\ and U2 are Holder continuous (in 
x) a priori. In fact this regularity follows from the above theorem! To understand 
why, and to see details about the derivation in the case O = R''^, we refer to [21] . 
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5.3. Some quasilinear equations. Consider equations of the form 

-ti-[a{x,Du)cr{x,Du)'^D'^u]- f{x,u,Du) + \u = Q in n, (5.4) 

where A > 0, f{x,r,p) continuous, increasing in r, a{x,p) = a{x,p)a{x,p)'^ , and 

\p-<l\ 



\a{x,p)-aiy,q)\<K(\x^y, , , ^ |^| ^ |^| 

\f{x,r,p) - f{y,r,q)\ < K[{1 + \p\ + \q\)\x - y\ + \p - q\]. 

In this case (HI) and (H3) hold in addition to an assumption similar to (H2). If 
we also assume (HO), (HBl), and (HB2), then existence and comparison for the 
boundary value problem (|5.4p and (|1.2p was proved in [5] . 

More general fully non-linear equations with "quasilinear" gradient dependence 
can also be considered. We omit this to get a shorter and clearer presentation. 
For the same reasons we also restrict ourselves to the case of Lipschitz continuous 
solutions and data, i.e. a = P = l,r] = e. 

In this case the quasilinear term in the equation gives rise to a term like 

\Wip)-aiq)\' 

in the proof of the comparison result (when a does not depend on x) . By (jA.15[) 
in Lemma |A.5[ 

\p^q\<K\p\A\q\\x~y\+K{e^ + B) 
when e is small enough, and hence by the assumptions on a, 

:\a{p) - a{q)f < —\x-yf + 



£2'^^^ '^'' - £2 1 ^1 e2l+|p|2 + |^|2 



<—\x-y\^ + Ke^ + —B^. (5.5) 



This computation motivates replacing assumption (Dl) by: 

(Dlq) There are 6i, 82 > 0, and Kf{K) > such that for any X > 0, 

F2{y,r,q,Y) - F,{x,r,p,X) < Kf{K)(^^ + <5i + 1522+5+1^2^ 

iorO<e<l,B>0,x,yeO.,reR,\r\<K,p,qe M^ and X, F e §^ satisfying 
\x-y\< Ke'^, \p-q\< K\p\ A |q|e2 + K{e'^ + B), \p\ + |<?| < K(l + e^ + B), and 

X \ ^ K f Id -Id\ , ,^, 2 , ON fid 

~y)^ 7^[-id wj+^^(^ +^)U Id 

The continuous dependence result now becomes: 

Theorem 5.3 (Quasilinear equations, Lipschitz solutions). Assume (HO), (HI), 
(H3), (HBl), and (HB2) hold for Hi,H2,Gi,G2, and ui,U2 G C"^^{n). Define 
Z/2 = (i.^ Vz/2)(j^i A 1/2). 

// (Dlq) and (D2) hold and ui and U2 are sub and supers olutions of the boundary 
value problems (j2.2p and p.3[) respectively, then there exist a constant C > 
(depending only on Kp, K , Kq) such that 



Amax(Mi-u2) < c(5i+82 + — + — 
Proof. The proof is similar to the Proof of Theorem 12. 21 with two exceptions: 
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(i) Assume 61,62,^,^ < C^^ where C is big enough (the general case follows 



since ui,U2 are bounded). Since s is chosen in terms of 61,62,^ 



Ml M2 



a 



suitable choice of C will ensure that e is small enough such that (JA.ISP of 
Lemma [A . 5 1 holds . This estimate is needed before one can apply (Dlq). 
(ii) At the end of the proof the following estimate will appear (remember rj = e) 

Xi[ui{x) - U2(x)) < K [rj +61 + —62^ ^ + ~( ) )' 

where the new final term in the right hand side of the inequality is a con- 
sequence of the -n-B^ term of (Dlq). Minimizing e like we did in Theorem 
12.21 then gives the result. □ 

As an example we consider an anisotropic quasilinear equation with capillary 
boundary condition. The type of non-linearity appearing here is similar to the 
non-linearity appearing in the mean curvature of graph equation. 



-tr 



rp / DU ® DU \ ^O 



+ Aw + /(x, u, Du) = in O, 



^-e{x){i + \Du\^Y'^ = Q in dn, 

on 

where 6 is Lipschitz continuous satisfying \9{x)\ < lo < \ and / satisfies the assump- 
tions mentioned above. Assume ui and U2 are Lipschitz solutions of this boundary 
value problem with different cri,(T2, ^1,^2 but with same / and A. Then we may 



apply 'i'heorem|5.3|with di = = /^i, ^2 = ^1 — ^20, and 


82 = sup 

pGK« 

to obtain 


r'^'V i+w) r''^ i+bp^ 

\.\ui-U2\ <c(|ai-a2| + -|^i-^2|o). 



Remark 5.1. Neither the above assumptions on a nor assumption (Dlq) is satisfied 
p-Laplacian equations. 

Appendix A. The test function 4>a 

A.l. Construction. The construction of the test function follows G. Barles [HH] 
tracking some quantities more precisely than he needs to. We give some more 
details (compared to [Tl[2]) whenever we feel this is helpful for the reader. 

First let d be a W^'°° extension of the (signed) distance function from some 
neighborhood of d^ to M^ (by (HO) the distance function is W^'°° near dVt). 
Furthermore we may and will choose d such that < d < 1 and \Dd{x)\ < 1 in 
Ct. We also extend the outward normal vector field n of dVt to all of M^ by setting 
n{x) := —Dd{x) for x e i7. An important consequence of the W^'°° regularity oi d 
and Taylor's Theorem, are the following inequalities: 

±[d{x) - d{y)] < ±iy - x) ■ n(^) + ^|i?'rf|o|x - y\\ (A.l) 

Next we extend Gi,i = 1,2, to a neighborhood V of dO. such that properties 
(HBl) and (HB2) still hold here (possibly with different constants K and v). Then 
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(HBl) and the intermediate value theorem that imphes that there exists unique 
sohitions Ci{x,p), of 

G,{x,p + a{x,p)n{x)) = 0, z = l,2, (A.2) 

for every x (in the neighborhood V) and p. Note that 

G^{x,p) - Gi{x,p + a{x,p)) = G^{x,p) = G^{x,Q) + {G^{x,p) ~ Gi{x,Q)) 

and 

Gi{x,p + Ci[y,q)n(x)) - Gi{x,p + Ci[x,p)n{x)) 
= G.i{x,p + Ci{y,q)n{x)) = Gj(x,p + C,(y,q)n(x)) - Gi{y,q + G^{y,q)n{y)), 

so by (HBl) and (HB2) 

ly^\C,{x,p)\<K{l + \p\) (A.3) 

and 

u,\C,{x,p) - C,{y,q)\ < K{1 + \p\ + \q\)\x-y\+K\p- q\ (A.4) 

for x,y G V, p,q £ M^. Now we extend Gi,i — 1,2, from V to R^ (in the x 
variable) such that (jA.4() and (jA.3[) are preserved (possibly with bigger K's). Next 
note that 

Gi{x,p + Gi{x,p)n{x)) - Gi{x,p + G2{x,p)n{x)) 
= -Gi{x , p + G2ix , p)n{x)) = G2ix,p + G2ix,p)n{x)) - Gi{x,p + C2{x,p)n{x)). 

Therefore, for x G 951, by assumptions {HBl) and (D2) we get 

i/i(Ci(x,p) - C2{x,p)) < Kg(hi + /i2 [\p\ + K{1 + IpI)]) . 
We have proved: 

Lemma A.l. Assume (HO), (HBl), and (HB2). 

(a) There exist unique functions Ci and C2 such that equation (IA.2P holds in a 
neighborhood of dfl, and the bounds (IA.3P and (jA.4|) hold for all x,p £ R^ . 

(b) If in addition (D2) holds, then there exists a constant Kc > such that for 
every p G R^ and x G dil, 

Gi{x,p) - C2{x,p) < -^(/^i + /i2(l + bl)). (A.5) 

Vl V 1^2 ^ ' 

To build smooth functions from C\ and C2 we will use the following sophisticated 
regularization due to G. Barles [2]: 

Definition A.l. Let a > 0. For a function C{x,p) on R^ x M^ we define 

where 

k = {a^ + {p-n(x)f)^ and r=(l + |pn^, 
and p G C°°(K.^) is a nonnegative function with total mass 1 and support in |a;| < 1. 
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Lemma A. 2. If C{x,p) satisfies (|A.4[) for x,p G M^ and < a < 1, then for any 

x,peR^, 



\Caix,p)\<Kr and \Ca{x,p) ~ Cix,p)\ < K{a + \p ■ n{x)\), (A.6) 

\D,Ca{x,p)\<Kr, \DpCaix,p)\<K, (A.7) 

^, \D^pCaix,p)\ < Kj, \DppCa{x,p)\ < ^, 



\D,.Ca{x,p)\ < K—, \D^pCa{x,p)\ < K-, \DppCa{x,p)\ < -. (A 



The proof follows from the classical properties of convolution and the regularity 
of C (|A.4p together with the choice of A and T. 

Now remember that d is a W^'°° extension of the distance function, and let 6*2,0 
be the smooth function obtained from Definition lA.ll with C = C2 . The full test 
function takes the following form: 

Definition A. 2 (The test function (f>a)- 

<i>a{x, y) = l|a; ~ j/p + ^ {d{x) - d{v)f - B{d{x) + d{y)) 
x + y 2{x~y) 

where A,B>Q are constants. 

A. 2. Properties. In the next 3 lemmas we the state main properties of the test 
function (f>a that we need in this paper. 

Lemma A. 3 (Positivity). Assume (HO), (HBl), (HB2), and let (j> be defined in 
Definition \A .2\ If A is big enough (not depending on e,a,B), then 

M^: y)>^\x- y\^ - Kos^ ~ B{d{x) + d{y)). 

Proof. By ((X6)) . \C2m{^,p)\ < C{1 + \p\) with C independent of a. So by Young's 
inequality one can take A big enough to insure that 

^1^ - yP + ^ (rf(^) - diy)f - C2A^, ^^^^){dix) - diy)) > -Koe^ 
which proves the Lemma. D 

Lemma A. 4 (Boundary conditions). Assume (HO), (HBl), (HB2), (D2), < 
£j '?! (J f^ 1; o.iT'd let (j)a be defined in Definition \A.2l Then for any x,y G Cl such that 
2^ ~ J/| < KirjE, there exists a K > only depending on the data and on Ki, such 
that if 

„2 , 2 , ^^ , K ( .. , .. "n 



B ^ K(n^ +e^ + a)-\ M1+M2- and A^K, 

1^1 V 1^2 ^ £'' 

then 

Gi{x,D,,4>a{x,y))>0 ifxedn, (A.9) 

G2{x,-Dyci>a{,x,y))<{) ifyedQ. (A.IO) 

Proof. We only prove (|A.9p . the proof of (jA.lOp is similar but easier due to the 
choice of C2,a in the test function. Note that d{x) — 0, d{y) ~ d{x) = \d{x) — d{y)\, 
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and remember that n = —Dd. We have 



Dx(l)a{x,y) 



\DxC2AX.p) + \DpC2,a{X,p)\ {d{y) - d(a;)) 



2A, 



+ C2.a{X,p) + B + ^{d{y)-d{x)) n{x) 
= p + Ri + [Ci{x,p)+R2]n{x), 
where p — -^\x — y\, X ^ ^{x + y), and 



(A.ll) 



2A 



R2 := hCi(x,p) + C2{x,p)] + [-C2{x,p) + C2,a{X,p)] +B+ -^{d{y) - d(a;)) 



According to (|A.6p . (|A.7p and (jA.ip . the second term in R2 is bounded below by 

Here we have also used that p — -^\x — y\ and |a: — y| < Kerj. 
According to (jA.5p . the first term in R2 is bounded from below by 



Kc 



(Mi+M2(l + bl) 



K f 7/ 



This means that i?2 > if B > K{ff + e'^ + a) + -j^ (fii + /xaf ) and A > f . By 
(|A.7|) of Lemma [A. 21 the regularity of d, \x ~ y\ < Ker], we also find that 



2 , ^2 , Kx) -d{y)\ 



\Ri\ < [Kil + \p\) + —K^ {d{y) - d{x)) < K(^r]^ + e^ + 

By (HJBl), (HJB2), (|A.11|) . (|A.2p . and the above estimates and choices of A, B 
Gi{x,Dx(j>a{x,y)) > Gi{x,p + Ci{x,p)) + i^iR2 - K\Ri\ 

2 , ^2 , \dix) -d{y)\ 



B-Ki,^+s^ + a)-^U+,2'-)+i2A-K) ^'^'^-/^y^^ 



K\t^' +e^ 



and the right hand side is strictly positive if 



viB = {l + vi)K{rf- +e^ + a) + — (^^1+^2-) and 2viA = {I + vi)K. 

v\\J V2\ eJ 

D 

Lemma A. 5 (Derivatives). Assume (HO), (HBl), (HB2), and let x,y ^ 0. and (pa 
be defined in Definition \A.S\ . Then 

\Dxd)a{x, y) + DyMx, y)\ < J^( '^~/' + e^) + 2B, (A.12) 

\D,Mx, y)\ + \DyMx, y)\ < K{e' + ^^ + ^^) + 25. (A.13) 

Furthermore if -^\x — y\'^ < C and A big enough (independent of a,e, B) then 



\'Dx^a{x,y)\,\Dy4>a{x,y)\ > 



2£2 



(1 ^s^B + K{1 + A)]) -K{e^ + B), 



(A. 14) 
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and if in addition e^ < [2B + K{1 + A)] then 

\D,Mx:y)+DyMx,y)\ < K\D.,M^:yMDyMx,y)\\^-y\+K{e^+B). (A.is) 

Finally, if \x — y\ < Krje and < 77, e, a < 1 then 

Proof. To simplify the computations, we make a change of variables and define $a ^ 

by 

'/'a(x,y)=$a,e(^,>^,^,r) 

where X = ^, Y^x~y, Z ^d{x)-d{y), T = d(x) +d(2/). We see that 

$a,s(^, y, Z, T) = ^ - C2,a(^, ^)Z + ^Z2 - ST, (A.17) 

e £ £ 

and straightforward computations lead to 

2Y 

z?y<i>,,,(x,y,z,r) = — - _i?^c2,.(x,— )z, ^^^^^^ 

2y 2v4 



and 



^x0a(a;, y) =]^Dx^a^e + Dy^a^e + Dz'i>a,eDd{x) + DT'i>a,eDd{x), (A.19) 



i?y(/.a(x, y) =-Dx<^a,e - DY<i>a,e + D z<^aA- Dd{y)) + DT<^a,eDd{y). (A.20) 



First note that estimate (JA.ISP easily follows from (|A.12[) and (jA.14[) . We start 
by proving estimate (|A.12p . Using the fact that n — —Dd, we see that 

D^(f>a{x, y) + Dy(f>a{x, y) = Dx^a,e + Dz<^a,e{n{y) - n{x)) + B{n{x) + n{y)). 

The first term on the right hand side can be estimated by (|A.7[) of Lemma IA.21 

\Dx^aA < i^(l + ^) \d{x) - d{y)l 
and the second term by (jA.6|) . 

C2,a(^,§^)|<i^(l + ^) and \Dz'^aA<K[l+^^) + '^\d{x)~d{y)\. 
Combining these estimates, using the regularity of d, and Y = x — y then leads to 

\D^cl>a{x, y) + Dyc^aix, y)\<[K + 2A\Dd\a) (l +^^)\x-y\ + 2\Dd\oB, 

and estimate (jA.12[) follows from Young's inequality. In a similar way, we can also 
prove C03l) . 
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We proceed to prove the lower bound ()A.14p and start by estimating \Dx4ia{x^ y)- 
Dy(j)a{x^y)\. Observe that 

2(x — v) 2 2Y 

D^cjyaix, y) - Dy^^, y) = 2[ ^ ^^ ^' - -DpC2AX, -^){d{x) - d{y))] 

2Y 2A 
+ {Dd{x) + Dd{y)) [ ~ C2,a(X, -^) + -^{d{x) - d(y))] - B{Dd{x) - Dd{y)) 

Now using a Taylor expansion and regularity of d, we see that 

{Dd[x) + Dd{y)) ■{x~y)< 2{d{x) - d{y)) + \\D''d\o\x - y\' , 

and after applying also (jA.6|) we get 

{Dx(l)a{x, y) - Dy(j>a{x, y)) ■ {x - y) 

It — 7;P 4 2Y 

^ ^-^ - ^DpC^AX, -72) -i^- vMx) - d{y)) 

9Y A A 

+ 2C2AX, ^){d{x) - d{y)) + —{d{x) - d{y)f 

- \D^d\,\x - y\^ (^\Dd\^\x -y\ + K{l + ^-^ 

-B{Dd{x)-Dd{y))-{x-y). 

Using Young's inequality as in Lemma IA.3I and taking A even bigger if necessary 
(but not depending on a, e, B), we have (e < 1) 

{D.x(l)a{x,y) - Dy(l)a{x,y)) ■ {x - y) 
> ^^ - B\Dd\o\x - y\' K{1 + A){1 + e^)^^ 



2- 



l^-^l (l-e^(B^K(X^A)^^^ 



and Cauchy-Schwarz inequality immediately yields 

\Dx<l>,{x,y)~Dy^,(x,y))\>^^^^{\-e'(B^K{\ + A)'^^ 

Now (IA.14P follows by combining the last inequality and (|A.12p . 

\Dx<i>a{x,y)\>\^^{\-e^{B + K + K{\ + A}^^^^-K[e' + B). 

Now we prove estimate (|A.16|) . A straightforward calculation using (JA.ISP yields 

2Y 

£^ 

2Y 
Dxz<i>aAX, Y, Z, T) = ~DxC2AX, ^), 

2 2Y 

DxY<^aAX,Y,Z,T) = --D^pC2AX,—)Z, 

e^ e 

2 / 2 \ 2 2Y 

DYY^aAX,Y,Z,T) = -Id^ (-j DppC2AX,—)Z, 

2 2Y 2A 

DYZ^aAX,Y,Z,T) = --I?pC2.a(A,— ), Dzz^aAX,Y,Z,T) = —Id, 
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We will estimate these terms using Lemma IA.2I and (|A.1|) . For example, one has 
\Dxx-^aAX,Y, Z,T)\ < K^[d{x) - d{y)\ < ^^(^'Ip ' «WI + \D^d\o\x ~ yj^), 

where p = ^- In our case A = [a^ + {p ■ n(A))^] 2 , and F = (1 + |pp)2 , and hence 

F^ 

.2t^2 ^ 



iDxx'^aAX, Y, Z,T)\< Kie^T^ + —\x - y 



a 

<K{e^ + \^){l+'-\^). (A.21) 

By carefully doing computations like above, one can prove that 

\Dyy<^,AX,Y,Z,T)\ < ^{l + it^)^ (A.22) 

\DYz'^a.e{X,Y,Z,T)l\Dzz<^aAX,Y,Z,T)\ < ^, (A.23) 

\Dxz<i>aAX,Y,Z,T)\<K{l+^-^^-^) (A.24) 

\DxY<^aAX,Y,Z,T)\ < K{\ + ^^)(1 + ^^^)- (A.25) 

Now we compute the matrix D'^(pa{x, y) from (IA.19[) and (|A.20|) : 

D^(l)aix, y) = Ml + M2 + M3 + MA + M5 + A/g + Mj + Mg, 

where 

^''^2 - T n2 d. n2 



4 Vi'ix^a.e ^ix*a 

l?iv$a.. 

2 



-i?iv$o 



Af -n2 d, f Dd{x)®Dd{x) ~Dd{x) (S Dd{y)\ 
^ - ^zz Va.e I _£)d(y) ^ i:,d{x) Dd{y) ® Dd{y) ) 



Mr. =1?^* 



2 



2Dd(a;) -^^(a;) - Dd{y)\ 



ZY-^a,e '■<y \ _ud{x) - Dd{y) 2Dd{y) J 

M -n^ <i> 6^1 ^^(^) \{Dd{x) ~ Dd{y))\ 



M7=Z?Z$a,5l n n2 



Mg = - S I „ „2 



D2d(a;) 

-D^d{y)J ' 

1)2^(0;) 

D^d{y)^ 



It can easily be seen that Mi (use (jA.22[) ). A/2 (use (IA.2ip ). and Afg can be bounded 
from above by a matrix of the form (IA.16P . Note that 

(A/3(C,^),(C,«)) = (^Xy*a,.(C-«),(C + «))<^|i?iy<i>a,enC-'^P + e'lC + «:P, 
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where 6 — r\^\ + ^rf- and hence by (|A.24|) Mj, is also bounded from above by 
(|A.16|) . Now we write 

D'^d{x) \ /O 



and handle the first part of M7 like we did with M3 . The second part can handled 
using the Ty'^'°°-regularity of d together with the first order estimates of <&a,e- We 
proceed with M4: 

(Af4(C, «), (C, '^)) = ^ (C • Dd{x) ~ K ■ Dd{y)Y 

<^4\Dd\l\C^ .\' + \D'd\Wm- + \n\'). 

The two remaining terms can be treated analogously using also (|A.22[) and (|A.24p . 
This ends the proof of the Lemma IA.5I D 
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